The absence of fragmentation in Bose-Einstein condensates 
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A Bose-Einstein condensate produced by a Hamiltonian which is rotationaUy or translationally 
symmetric is fragmented as a direct result of these symmetries. A corresponding mean-field unfrag- 
mented state, with an identical energy to leading order in the number of particles, can generally be 
constructed. As a consequence, vanishingly weak symmetry-breaking perturbations destabilize the 
fragmented state, which would thus be extremely difficult to realize experimentally, and lead to an 
unfragmented condensate. 
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Introduction. One of the most fundamental issues in 
the problem of Bose-Einstein condensation is that of frag- 
mentation. Years after the prediction of this phase tran- 
sition by Bose and Einstein, Penrose and Onsager [ij 
provided a rigorous criterion for the existence of a Bose- 
Einstein condensate. Starting from the TV-particle wave 
function of a system, one determines the eigenvalues of 
the (Hermitian) one-body density matrix, 

pW-ri) = 



TV 



dv2... drAr**(r'i,r2, . . . r7v)*(ri, ra, . . . ,rAr). (1) 



If at least one eigenvalue is of order TV, the system is Bose- 
Einstein condensed, otherwise it is not. A single eigen- 
value of order N indicates simple condensation; when 
more than one of the eigenvalues are of order N , the 
condensate is said to be fragmented [41. 

Years after this definition was introduced, Nozieres and 
Saint James [3| argued that, in Hartree-Fock approxi- 
mation, the Fock term makes it energetically favorable 
for the system to fragment if the effective interaction 
between the bosons is attractive. The question of con- 
densate fragmentation remained academic, since homoge- 
neous systems with an effective attractive interaction are 
unstable against collapse. Modern techniques for dealing 
with trapped cold atoms have rekindled interest in this 
question since these gases can be metastable if the effec- 
tive interaction is attractive, and it may thus be possible 
to realize a fragmented state. This issue has been the 
subject of a number of studies over the last decade, see 
e.g. , Refs. I, i i 0, i, S [M [O, E d • 

In the present study we consider the problem of frag- 
mentation for an effective attractive interaction between 
the atoms. As two characteristic examples we consider 
bosonic atoms at zero temperature that are confined in 
toroidal and in harmonic traps. We will argue that the 
issue of fragmentation is subtle. For any fragmented 
state, it is always possible to construct a mean-field, 
non-fragmented product state that has the same energy 
to leading and subleading order in N as the fragmented 
state. In these representative problems and in any case 
where the Hamiltonian is rotationaUy or translationally 



invariant, the single-particle density matrix is diagonal, 
and its eigenvalues are the occupation numbers of the 
corresponding single-particle states. The system is frag- 
mented merely as a consequence of the symmetry of the 
Hamiltonian. Due to the low excitation energies of states 
of well-defined (angular) momentum [i.e., characteristi- 
cally ©(l/TV)], such states are fragile and virtually im- 
possible to realize in practice. We shall show that small 
symmetry-breaking terms in the Hamiltonian that vanish 
in the limit TV — > oo can materially alter wave functions 
and can reconstruct an unfragmented condensate. 

Bosons in a toroidal trap with an effective attractive 
interaction. Consider a tight toroidal trap where the in- 
teraction energy between the atoms is much smaller than 
the excitation energy in the transverse direction. The 
Hamiltonian of this system reduces to [ij, llH 



H 



N 

E 



80? 



1 



N 



v{e^)--\u\ J2 m~o,), (2) 



i#i=l 



where 9 is the angle in cylindrical polar coordinates, V{0) 
is the potential that acts along the toroidal trap, and U 
describes the coupling and is proportional to the scatter- 
ing length for elastic atom-atom collisions. The dimcn- 
sionless parameter 7 = —(TV — l)|[/|/(27r) gives the ratio 
between the interaction and the kinetic energy. As shown 
in Refs. [IJ, [l^, the gas becomes localized and forms a 
bright solitary wave when 7 < 7j, = —1/2. 

It is convenient to expand the eigenfunctions of this 
Hamiltonian in the basis of states exp{iq9) / ^/2Tr , with 
q = 0,±1,±2..., which are eigenstates of the angular 
momentum operator, L — —id/dO, with eigenvalues q. 
For clarity, we truncate the space and consider only single 
particle states with q = —1, 0, and 1. We also use the 
notation for the corresponding Fock states 



|m) = |(-l)'",0 



N-2m-L 



,(+1)' 



(3) 



where to, N —2m— L and m+L are the occupancies of the 
single-particle states with q — —1,0, and 1, respectively. 
Clearly, the above states are eigenfunctions of the angular 
momentum operator L and of the number operator TV. 



The Hamiltonian of Eq. ^ commutes with N. If ¥(9) 
is zero or constant, it also commutes with L. For a spe- 
cific N and L, the eigenstates \iPl) of this Hamihonian 
can be expanded in the basis of the states |m), 



li'L 



E'^^K-ir'O 



N-2m-L 



A+i) 



'm-\-L\ 



(4) 



The eigenvalue equation H\iPl) — E{L)\iPl), where E{L) 
is the eigenenergy, has the form 16] 



^m^7n^rn ' ^ni.m- 



1< 



H, 



i< 



7n,m4'l"'rn+l 



EiL)dt{5) 



where _ff„ 



(nlHlm) are the matrix elements of the 



Hamiltonian between the states \n) and \m). In the spe- 
cific truncated space, only the matrix elements Hm^m and 
Hm,m±i are nonzero. In the limit N -^ cx), numerical so- 
lution to this equation yields 






0.488 +^^+0(l/iV). (6) 



For large N, it is convenient to regard m as a con- 
tinuous variable. Assuming that d^ is a differentiable 
function of m up to any order, 






it Oin^rii 



(l/2)af„d^ 



(7) 



and the matrix eigenvalue equation can be written as a 
familiar harmonic oscillator problem. 



We now consider the efi^ect of a very weak inhomoge- 
neous potential, V{9), on the eigenvalues of the single- 
particle density matrix. (This point has also been ad- 
dressed in Ref. i9|.) As we will see, the state is no longer 
fragmented in the presence of such an inhomogeneity. We 
choose V{9) ex 6 cos 9, with S ^ 1. This potential con- 
nects single particle states with Aq = ±1, 

V{9) = d{aoai + aiaj + aoaLi + a-iaj), (10) 

with Gq and aj the usual annihilation and creation op- 
erators of particles with angular momentum q. For suf- 
ficiently weak V{9), the resulting eigenstates jV') of the 
Hamiltonian, Eq. ([2|), can be expressed as a linear super- 
position of the states |?/'l)j 



E^-K-ir'O 



N-2r, 



\Tn+L\ 



(11) 



m.L 



Given the form of V{6) , it can connect the state 

i(-ir 



0^ 


-2m-L^ (-1-1)™+^) to the following four states 


ll> 


= \i-i)"'-\o^-^"'-^+\{+ir+^), 


|2) 


= |(-i)"+i,o^-'™-^-\(+i)"+^). 


|3) 


= |(-i)™,o^-'™-^-\(+i)"+^+i). 


|4) 


= |(-l)",0^-'™-^+\(+l)™+^-i). (12) 



Thus, Eq. ^ assumes the form, 



1 
2/i 



d^ d^ 



El^ + ^{m-m^)^ 



E{L)di. (8) 



The solution of this equation, with the boundary con- 
dition that d^ vanishes as m ^ zLoo, has the expected 
Gaussian form 



(9) 



d^ ex exp [-^yXfi{m - mf^Y /2\ 



-1 we find that TOq 



For the specific value of 7 
0.1429iV - L/2, A = 14/iV, and l/fi = 0.20577V. The 
root mean square deviation of m from tjiq scales like yN. 
The number of states participating actively in Eq. ^ is 
0{y/N), and the approximation of Eq. ([8|) is justified in 
the large N hmit. Similarly, we find Eq « -9iV/14 - 
0.357 -I- {3/2)L'^/N. Our approximation to the energy, 
E{L) — Eq + {\/2)^/JJJi, thus agrees with the exact 
result of Eq. ^ for all terms shown. The error in energies 
obtained with this smooth approximation is 0{1/N) for 
all \L\ < O(ViV) to be considered below. 

The one-body density matrix corresponding to the 
state \4>l) described by Eqs. ^ and ^ is diagonal in 
the basis of angular momentum eigenstates. This is a 
direct consequence of the fact that the states \iPl) are 
eigenfunctions of the total angular momentum. The di- 
agonal elements are the occupancies of the three single- 
particle states with q = —1,0 and 1. Since the average 
value of m = mo is of order A'' (but not equal to N), 
all three eigenvalues are of order A^. The condensate is 
fragmented. 



^m,m,^rn 



H, 



-5^JmQ{N -2mQ){d, 



m , m - 
L+1 



-id 



-1 



H, 



m ,rn-\-l^r-n + l 



-lL-1 



Edi 



jL+1 



(13) 



Here, ^JmQ{N — 2too) approximates the value of this 
prefactor with its value at the minimmn, i.e., mo = 
771^=0 = 0.143153A. 

Converting this matrix eigenvalue equation into a dif- 
ferential equation as above, we find 



1 
2^ 



d^ d^ 



A 



L\2 



£^0 +77("i-0 



+ 5^/m^{N^^2mo){'id^ + 2^LLd^ + 



d„ 



2dLdmdi:,) = Edt 



(14) 



The solution of this differential equation has the form 



d„ ex exp [— ai(m + L/2 — tuqY — 02^^ 



(15) 



with ai and 02 positive, since this function vanishes as 
m and L tend to infinity in any direction. Numerical 
calculations with a symmetry-breaking term V{9) of the 
form (5 = -(l/100)(100/Af)^-i5 verify that the coeflicients 
d^ are indeed Gaussian distributed as a function of m 
and L. 

Direct calculation reveals that the one-body density 
matrix now has only one eigenvalue of order A^. [The 
next-largest eigenvalue is 0(A^°-^™).] This result is read- 
ily understood as a consequence of the Gaussian support 



of the d^. If the symmetry-breaking potential, V{9), is 
sufhciently strong that the root mean square variations 
in TO and L grow with A'" but small enough that Am/{m) 
and AL/{L) vanish as iV ^ cx3, the matrix elements of 
the one-body density matrix, pij = {alaj), are ^Jn^rTj, 
with Ui the occupation number of single particle state i. 
The one-body density matrix is thus rank one separable. 



It has one non-zero eigenvalue of ^ 



iV; the ele- 



ments of the corresponding eigenvector are proportional 
to ^JfTj. All other eigenvalues are zero, and the conden- 



sate is unfragmented. 

The symmetry-breaking potential, y{Q\ must be suf- 
ficiently strong if it is to yield the desired mixing of the 
states |'0l)- According to Eq. ©, the states with L 7^ 
are separated from the L = ground state by a term 
that scales as I? jN . Since the contribution of V{Q) to 
the energy is of order A^(5, b must vanish less rapidly than 
1/iV^. In addition, low- lying excited states (for each L) 
are separated by an energy of ©(iV") from the lowest en- 
ergy state. Validity of the truncation to the states |-0l) 
requires that b vanishes more rapidly than 1/iV. Clearly, 
the second condition is dictated by approximations made 
in this calculation and not with the absence of fragmenta- 
tion. In short, even symmetry-breaking potentials which 
vanish in the large N limit are sufficient to ensure that 
the condensate is not fragmented. 

Bosons in a harmonic trap with an effective attrac- 
tive interaction. We now turn to two additional systems 
which have been considered as examples of condensate 
fragmentation. First, consider rotating bosonic atoms 
confined in a two-dimensional harmonic trap and subject 
to the Hamiltonian 4]. In cylindrical polar coordinates. 



^ 1 



V 



2 , 1 2 1 



i=l 



N 



r,). (16) 



Here 77, which describes the atom-atom interaction, is 
proportional to the s-wave scattering length. If the cou- 
pling is weak, jryj <C 1, this Hamiltonian can be trun- 
cated to include only states in the lowest Landau level 
with zero radial nodes and m quanta of angular momen- 
tum, <^m = 2;™e~'^' '^/vttto!, with z = x -\- iy, where 
X and y are the Cartesian coordinates. In this case, as 
shown by Wilkin, Gunn and Smith 3 and by Mottelson 
[17| , the interaction energy of the lowest energy state for 
any given angular momentum is the same as that of the 
non-rotating system. 



E{L) = -^-^N{N 



\M'd^P^~^^N{N-l){n) 



The full energy of these states contains an additional con- 
tribution of |L| -I- 1 from the confining potential. The cor- 
responding exact many-body eigenstate describes a cen- 
ter of mass excitation with 



N 



^.^.(zi, Z2, . . . , zn) = NlZ"^ V exp(-|z,| V2). (18) 



=1 



I/Vtt^N'^LI and Z is the center of mass 



'7 V^^ 



Here, A^i = 
coordinate, i.e. 

The eigenvalues of the single-particle density matrix 
are p,„ = {N - 1)^-"'L\/[N^-\L - to)!to!] Q. Due 
to the axial symmetry of the Hamiltonian, this density 
matrix is diagonal, and its eigenvalues are simply the oc- 
cupation numbers, iV|c„ip, of the single-particle states. 
The energy is minimized when all Cm have the same 
phase, which can be taken as positive without loss of 
generality. In the limit of infinite N and L with / ~ L/N 
finite, we see that [18] |cmp(/) = f"exp(-/)/TO!. Ac- 
cording to the usual criterion, this is a fragmented state. 

It is possible, however, to construct a mean-field, prod- 
uct wave function which has the same interaction energy 
and which is necessarily unfragmented. Consider the sim- 
ple form 



N 



*W(^1>^2,. 



= TT V" C„i(/)m(zi), 



i—1 m—0 



(19) 



with the coefficients Cm — ^yi"'^/mlexp{~l/2). This state 
is normalized, and the expectation value of the angular 
momentum per particle is I = L/N . The interaction 
energy of this state can be calculated analytically, 

<F = -^N{N-l) j \Y, c„r4>m{z)\^dxdy 



in 



N{N-l)e~^^J2 



m=0 
1 / / \ "^ ^^ 

i I I \ \-^ I m \ I m 



E 



to! V2/ '^ \ k ) \ J 

m=0 ^ ^ k,j=Q ^ ' ^ 



W _._ .. _9, -sr^ 1 f I 



^"("-^-^'EiU ^ 



47r 



47r 



iV(iV-l), 



(20) 



which is identical to the energy given by Eq. (J17p . The 
overlap between the states ^^^ ^.nd ^mf '^a.n also be cal- 
culated analytically as 



(*c';I*Lf) 



Ne) 



\/tt'^N^LI V2^ 



(21) 



This overlap vanishes in the thermodynamic limit L -^ 
00. This comes as no surprise, since ^^^ describes a state 
that is spread uniformly around the center of the trap, 
while ^I'mf describes precisely the non-rotating "clump" 
of matter displaced from the center of the trap and ro- 
tating around it. 

The existence of a mean field state with an energy 
close to the exact eigenvalue is a relatively general conse- 
quence of rotational or translational invariance. Given a 
Hamiltonian which is axially or translationally invariant, 
the one-body density matrix is diagonal with eigenvalues 
equal to the occupation numbers of the corresponding 
single particle states. From these occupancies, it is possi- 
ble to construct a mean-field wave function with the same 



energy as the exact solution to leading and often sublead- 
ing order in the number of particles, N . The same conclu- 
sion applies to the argument of Nozieres and Saint James: 
The exchange interaction does not necessarily favor frag- 
mentation since there exists a non-fragmented mean-field 
product state with the same energy in the iV — » oo limit. 
As seen in the case of toroidal confinement, the question 
of whether the system is better described by a wave func- 
tion which is an eigenfunction of the total momentum or 
angular momentum (and thus fragmented) or is better 
described by a mean field wave function (and thus not 
fragmented) must depend on the response of the system 
to vanishingly small symmetry-breaking terms. 

This issue can be investigated with arguments and con- 
clusions identical to those above. A one-body symmetry- 
breaking term, V{z) = 6(z + z*) is introduced, where 
again z = x + iy. The basis of states is truncated to in- 
clude only the lowest-energy states, *^^. Evidently, V{z) 
can only connect the state L with the states L ± 1. The 
corresponding matrix elements are, e.g., 



{^^+'\V\^^J = 6^N{L + 1) for i > 0. 



(22) 



If 6 vanishes with increasing N, this truncation of states 
is legitimate. If it vanishes more slowly than 1/\/N, there 
will be significant mixing of the states 'i'^^. As in the 
case of toroidal confinement, the wave function will have 
localized (i.e., Gaussian) support in the space of single 
particle states. Precisely as before, only one eigenvalue 
of the one-body density matrix is of order N, and the 
condensate is not fragmented. 

Identical arguments can be applied to the related but 
simpler two-state model of Ref. [3|- There, the conden- 
sate is fragmented due to a "parity" symmetry. This is 
reflected in the fact that eigenstates contain, e.g., only 
an even (or odd) number of particles in one of the states. 
The energy difference between the lowest-energy even 
and odd states vanishes exponentially with N. Once 
again, a vanishingly small one-body symmetry-breaking 
term, proportional to (ajai + a[ao), is sufficient to recon- 
struct an unfragmented condensate. These two examples 



will be described in greater detail elsewhere. Finally, sim- 
ilar arguments apply to the studies of Refs. [l3,[i3l> where 
a state fragmented by some symmetry of the Hamiltonian 
can be restored to a simple unfragmented condensate by 
very weak symmetry-breaking perturbations. 



General conclusions. For many systems of a large 
but finite number of bosons with attractive interactions, 
mean field theory provides a good description of the 
ground state energy and leads to the unambiguous pre- 
diction of an unfragmented condensate. The imposition 
of general constraints, such as conserved total momentum 
or angular momentum, characteristically produces min- 
imal changes in the ground state energy and frequently 
indicates condensate fragmentation. This apparent con- 
tradiction has led some authors to suggest modified cri- 
teria for condensate fragmentation. We have offered an 
alternate resolution. The various examples considered 
here all suggest that the small excitation energies of ex- 
cited states in these systems can render them sensitive 
to vanishingly small symmetry-breaking perturbations. 
The resulting localized (i.e., often Gaussian) support 
of the wave function then leads to a one-body density 
matrix, approximately given as pij — ^/uifTJ, which is 
rank one separable with one eigenvalue of 0{N) in the 
N —> oo limit. The unfragmented condensate, decon- 
structed by rigorous symmetries, can be reconstructed 
by small symmetry-breaking perturbations. Such per- 
turbations can be difhcult to eliminate experimentally. 
(Such antagonism between the mean field approxima- 
tion and symmetries is well known. For example, insis- 
tence on maintaining translational invariance in fermion 
systems leads inevitably to a trivial Hartree-Fock wave 
function of plane wave states and a poor description of 
both the ground state energy and wave function.) While 
the present results in no sense rule out the possibility of 
condensate fragmentation in systems of bosonic atoms, 
they do suggest that it is important to demonstrate that 
theoretical indicators of fragmentation are robust with 
respect to small symmetry-breaking perturbations. 
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